publication ^ artkle SUbmitted '^journal for 

publication' rhis ° fchanges that ma y b e made before formal 
that it w M nm h PreP ',' nt 15 made avai,ab]e the understanding 
Author C ' ted " repr ° dUCed Without the Permission of 1 


A/ 7~ 


Gaussian beam intensity flattener 


DeVon W. Griffin 

NASA Lewis Research Center 
Microgravity Fluid Physics Branch 
M/S 100-3 

21000 Brookpark Road 
Cleveland, Ohio 44135 
E-mail: devon.griffin@lerc.nasa.gov 

Subject terms: lens design; Gaussian beams. 

Paper CM-2127 received Dec. 24, 1997; revised manuscript 
received March 14, 1998; accepted for publication March 14, 1998. 
[S0091 -3286(98)03707-6] 


1 Introduction 

A recent article 1 outlined the procedure for using a gradient 
index lens system to change a Gaussian intensity distribu- 
tion to one that is uniform without any loss of flux to the 
Me 1 limit. While design of this system was straightforward, 
practical implementation has been delayed due to fabrica- 
tion difficulties. 2 Even though the pioneering work on such 
intensity transformations was performed nearly 30 years 
ago, 3 implementation difficulties have been the genesis of 
discussions involving many different approaches such as 
aspheres, 4 holograms, 5 and aberrated spherical optics. 6 





COMMUNICATIONS 


Shafer 6 discussed the use of spherical elements to generate 
a uniform beam to the Me diameter. The goal of this inves- 
tigation was to use commercial elements and extend the 
correction to a 1/e 2 diameter of 3 mm over long propaga- 
tion distances. 


2 Design Procedure 

The starting point for this design consisted of using two 
doublets, as suggested by Shafer. 6 Rather than use spherical 
aberration as a tool to redistribute the energy, as he and 
Frieden 3 did, ray locations were used to represent energy 
distribution as discussed by other workers. 1,4,5 In those 
cases, however, use of either custom gradient or aspheric 
optics assured perfect correction. In this investigation, nei- 
ther option was used. While the net result for one of the 
designs presented here was excellent to the Me point and 
acceptable to the Me 2 points, the results do suffer from the 
inability to independently address each zone, a capability 
that is only available with holograms, aspheres, or custom 
GRIN devices. 

The rayset used for optimization was chosen by first 
requiring that rays represent the boundary of a region of 
equal energy in the Gaussian beam. This stipulation re- 
sulted in a heavier sampling of the beam near the optical 
axis where the irradiance was highest. It was accomplished 
by calculating the volume under a unity amplitude two- 
dimensional Gaussian and dividing by 20, the desired num- 
ber of sampling locations in the merit function. A numeri- 
cal integration was then iteratively performed to identify 
the boundaries of these regions of equal irradiance. 

The targets in the merit function were the transformed 
locations of these rays in the output beam. The transforma- 
tion assumed conservation of energy between the input 
Gaussian and output uniform beams. Mathematically, the 
result is 5 


R = 


r z 

2^ [l-exp(-2r 2 /ro)] 


1/2 


( 1 ) 


where r is the radius of the input ray, R the position of the 
corresponding output ray, r 0 is the Me 2 radius of the 
Gaussian beam, and a is the uniform irradiance for correc- 
tion to that radius. 

Shafer 6 proposed four different design forms, but only 
extensively investigated the (+— ,— +) case to the Me ra- 
dius. When this form was used with the merit function 
previously described and correction extended to the \le 2 
radius, the beam remained approximately collimated be- 
tween the doublets while performing the energy redistribu- 
tion to what essentially was a first approximation. Though 
uniform illumination at a surface 50 mm following the last 
lens was achieved, that profile did not persist with increas- 
ing propagation distance. To increase the performance, 
both doublets were changed to triplets and the optimization 
continued. Once the change from the doublets had been 
made, the light was no longer approximately collimated in 
the center section. Rather, it formed a focus between the 
third and fourth and fourth and fifth elements. A split of the 
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Fig. 1 Layout of focusing system. 



Normalized Pupil 

Fig. 3 Irradiance distribution of both systems. 


first and second elements of the second triplet was required 
to obtain optimal performance. Figure 1 is a layout of this 
(+ + -,-) 1 1-) focusing system. 

An alternate design was found by using a negative lens 
as the first element with a catalog asphere following it in 
the first group. This avoided the real focus, which would be 
undesirable for high-power lasers. In this case the second 
element of the first triplet was eliminated since it did not 
contribute to improvement of the system. Additionally, the 
first element of the second group was split, yielding a de- 
sign of ( — +,H 1-—). As mentioned, the first element of 

the second group is aspheric, but it is a catalog element 
from Geltech. This system, hereafter referred to as the non- 
focusing system, is depicted in Figure 2. 

3 Final Design and System Performance 

In these designs, critical parameters were glass thicknesses 
and the distance between elements in the first doublet. Pref- 
erence on the former was always for thinner elements and 
the position of the latter was the most sensitive parameter 
in the design of the system which produced a focus. 



Fig. 2 Layout of nonfocusing system. 


The performance of the two design types is plotted in 
Figure 3. The desired input Gaussian and output uniform 
irradiance profiles are plotted for reference. Data points de- 
noted by darkened squares describe the performance of the 
nonfocusing system whereas the empty circles plot the ir- 
radiance for the focusing system. Both of these irradiance 
distributions were calculated 1 m from the rear vertex of 
the last element. As shown in this figure, the focusing sys- 
tem performs much better than the one that did not come to 
a real focus. The nonfocusing system may be viewed as a 
local minima in the search for the focusing system since 
great care was required to prevent the nonfocusing system 
from producing a focus. 

Figure 4 repeats the irradiance plot of Figure 3 for the 
focusing system but changes the range of the vertical axis 
to better display the irregularities. Additionally, empty tri- 
angles denote the irradiance of this system 2 m from the 
rear vertex of the last element. As illustrated, the beam 
changes very little over this propagation distance. Numeri- 
cally, the maximum error at both distances is 1.7% within 
the Me radius and a maximum of 2.6% at the edge of the 
field for the 1 m distance. That error grows to 4% at 2 m. 



Fig. 4 Irradiance distribution of focusing system at 1 and 2 m. 
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Table 1 Prescription for (+ + -,+- + -+) lens. Table 2 Prescription for (-+,+- + -) lens. 


Surface 

Number 

Radius 

Thickness 

Commercial Part 
Number 

1 

11.142 

2.9 

MG LPX 025 

2 


10.8438 


3 

5.19006 

1.9 

OS 011-0230 

4 


1.97153 


5 

-1.55 

1.0 

SH 314000 

6 

1.55 

189 

SH 314000 

7 

1 .42999 

1.5 

JMLCPX 10035 

8 


6.4156 


9 


0.5 


11 

1.55 

5.83376 

SH 314000 

12 

2.58399 

2.5 

SH 312007 

13 


0.67332 


14 


0.5 


15 

1.55 

2.76431 

SH 314000 

16 

-14.3 

3 

MF LFF088 

17 


2.38 

ES 45147 

18 

-4.65 




While the design could have retained use of custom ra- 
dii, catalog elements can greatly reduce the cost while de- 
livering adequate performance, as described elsewhere. 7,8 

Table 1 lists the prescription of the (+ + -,-! 1 1-) lens 

along with the associated manufacturer’s part numbers. 
Similar information is provided in Table 2 for the ( — +, 

H 1 — ) design. In these tables MG refers to Melles Griot, 

S/H to Spindler and Hoyer, NRC to the Newport Research 
Corporation, GE to Geltech, ES to Edmund Scientific, and 
OS to OptoSigma. All distances given are in millimeters 
and the designs assume a wavelength of 0.6328 gm. 

4 Conclusion 

The design forms for lens systems that will transform an 
incident Gaussian beam to an output beam with uniform 


Surface 

Number 

Radius 

Thickness 

Commercial Part 
Number 

1 

-12.92 

2.5 

NRC KPC 010 

2 


2 

MG LPK 004 

3 

7.77998 

0.6997'”“' 


4 

6.6276 

5 

CO 350220 

5 


200 


6 

16.962 

3.9 

MG LPX 057 

7 


2 


8 


2 

MG LPK 002 

9 

5.18699 

1 .9851 1 


10 


2 

SH 312201 

11 

-5.464 

3.93945 


12 

6.743 

2 

MG LPX 008 

13 


0.5 


14 

3.101 


SH 314001 


intensity over long propagation distances have been pre- 
sented. The quality of the correction at 1/e 2 points was 
presented. Prescriptions for implementation using stock 
lenses were provided. 
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